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We consider an s-wave superconductor in the vicinity of a second-order ferromagnetic (FM) or 
spin-density- wave (SDW) quantum critical point (QCP), where the superconductivity and mag- 
netism arise from separate mechanisms. The quantum critical spin fluctuations reduce the super- 
conducting Tc- Near a FM QCP, we find that Tc falls to zero as l/ln(l/K) in 3D and as k in 2D, 
where k ~ \ J — Jc^ is the inverse correlation length of the spin fluctuations, and measures the 
distance \ J — Jc\ from the quantum critical point. SDW quantum critical fluctuations, on the other 
hand, suppress Tc to zero as ^/k, in 2D, and suppress Tc only to a finite value in 3D, producing a 
cusp of the form const -|- | J — Jc]". 

PACS numbers: 



I. INTRODUCTION 

There has been much interest recently in the interplay 
of superconductivity and magnetic phenomena. Many 
different scenarios arise depending on the supercon- 
ducting pairing symmetry (singlet /triplet, s / p / d-wave)^ 
and the type of magnetism (static/fluctuating ferromag- 
netism/antiferromagnetism/spin density waves). Much 
attention has been given to magnetically-mediated un- 
conventional superconductivity, such as models of d-wave 
singlet superconductivity driven by antiferromagnetic 
spin fluctuations (possibly applicable to the cuprates and 
heavy- fcrmion supcrconductorsi) , and p- wave-like triplet 
pairing in He'^— i^, Sr2Ru04^ and UGe^^i^ driven by fer- 
romagnetic spin fluctuations. 

However, there are of course systems where conven- 
tional s-wave pairing coexists with magnetism. While 
the coexistence of antiferromagnetic or SDW order with 
superconductivity is not uncommon, the coexistence of 
FM with s-wave superconductivity is rare. In ErRh4B4i 
and HoMogSg^i^, magnetism and superconductivity arise 
from independent mechanisms (RKKY coupling between 
/-electrons and phonon exchange respectively). It is 
now well accepted that static magnetic order has a pair- 
breaking effect on BCS singlet s-wave superconductiv- 
ity (as does any perturbation which breaks time-reversal 
symmetry), often resulting in the complete destruction 
of superconductivity. The mean-field theory of ferro- 
magnetic superconductors was studied by Gor'kov and 
Rusinovi^, although a more general analysis suggests 
a Larkin-Ovchinnikov-Fulde-FerrcU state with a non- 
uniform order paramete r^^i^^ . A mean-field theory based 
on a simple 'nesting-fraction' model (similar to that of 
Bilbro and McMillani^ for the coexistence of charge- 
density-wave order and superconductivity) shows that 
coexistence of spin-density- wave (SDW) order and super- 
conducting order may also be possible in a limited range 



of parametersi^. 

One expects that magnetic fluctuations should also 
have a pair-breaking effect on s-wave superconductivity, 
especially for a system close to a second-order magnetic 
QCP, where the amplitude of the fluctuations becomes 
very large. Berk and Schrieffer— showed numerically 
that in the presence of ferromagnetic fluctuations a fi- 
nite electron-phonon coupling Vph is required to produce 
superconductivity, and in Pd, for example, this critical 
coupling is so large that superconductivity does not oc- 
cur at all. We study this problem in more detail, obtain- 
ing analytic results for the dependence of Tc as the QCP 
is approached. 

We assume that the pairing interaction is due to 
phonon exchange, and can be represented by a BCS in- 
teraction, whereas the magnetism is driven by a coupling 
J between itinerant electrons or localized spins such that 
when J exceeds a critical coupling at T = the sys- 
tem undergoes a second-order phase transition from a 
disordered state to a magnetically ordered state. The 
resulting interaction between electron spins (as a func- 
tion of Matsubara frequency iujn) is assumed to be of the 
'Hertzian' form 



Jkq^ 



aq 



or 



(1) 



near a FM(SDW) QCP. Here is the diverging corre- 
lation length. This approach is similar to the fluctuation- 
exchange (FLEX) approach used by Monthousi^ to study 
spin-fluctuation-mediated d-wave-like superconductivity, 
and to that of Wang et ali^ and Roussev and Milli&i^, 
who studied a model of p-wave pairing where the ferro- 
magnetic fluctuations suppressed Tc to a finite value at 
the QCP. Also, Li et ali^ have performed related stud- 
ies of the electromagnetic properties of superconductors 
near ferromagnetic instabilities. 
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As the QCP is approached, the energy scale of the spin 
fluctuations goes to zero ('critical slowing down') and 
their amplitude goes up. We show that in the vicinity 
of the quantum critical point, pair- breaking by slow spin 
fluctuations can be mapped onto the Abrikosov-Gor'kov 
problem of pair-breaking by static magnetic impurities. 
We thus predict that Tc may or may not be suppressed 
all the way to zero at the QCP, depending upon the di- 
mensionality and the type of magnetism: 



2D FM Tc 

2D SDW Te 

3D FM Tc 

3D SDW Tc 



K^lJ-Jcl" 

\J-JcY'^ 

l/ln(l/K) - l/ln(l/|J- 
const + I J - JcT 



where v is the correlation length exponent. We expect 
these forms to apply on both sides of the transition. 



II. CHARGE-CHARGE INTERACTION 

Starting with an electron-phonon Hamiltonian, inte- 
grating out the phonons in the harmonic approximation, 
and making certain approximations for the structure of 
the electron-phonon coupling in momentum space leads 
to the following action for phonon-mediated supercon- 
ductivity: 

S[ll>,-4>\ =T'^ {iSn - Ck)V'«ka^/'«kQ 

IK~t',^-^'Pt^Pt'^' (2) 

where e„ = 2ir{n+^)T are fermionic Matsubara frequen- 
cies, a = ±i is a spin index, is the electron dispersion 

relation, t runs from to /? = l/T, prx = X)a ^V'rxaV'rxa 
is the electron density operator, and is an effective 
interaction in the 'charge' channel that is second or- 
der in the electron-phonon coupling g and flrst order in 
the phonon propagator V. If the superconductivity is 
isotropic and the Fermi surface is spherical, or if the sys- 
tem is dirty enough that the order parameter A is effec- 
tively averaged over the Fermi surface, then for the pur- 
pose of calculating superconducting properties V^^ (the 
Fourier transform of V^^) '^^'^ be replaced by an effective 
interaction in the Eliashberg equations. is usually 
referred to in the literature as "a^F" , and is obtained by 
averaging the interaction V^^ over pairs of momenta on 
the Fermi surface^^: 



kp 



'5(ek) <5(ep) V, 



nkp 




'5(ek) 



(3) 



HI. EFFECT OF FM FLUCTUATIONS ON n 

A. Spin-spin interaction V,tq near a FM-QCP 

In the random phase approximation (RPA), the spin 
correlation function in an itinerant ferromagnet has the 
following form at small n and q: 



-v/O 2 



aq 



(4) 



where k, the inverse spin-fluctuation correlation length in 
the interacting system, is reduced from its bare value kq 
as K = kq (1 ~ J/ JcT = 1/^ with v = 1/2, so that x un- 
dergoes Stoner enhancement, k indicates the 'distance' 
from the ferromagnetic quantum critical point (QCP); a 
second-order quantum phase transition occurs at J — Jc^ 
when K = 0. Eq. ([4]) can also be obtained by integrating 
out the electrons to obtain a Landau-Ginzburg- Wilson 
effective action for a spin field. 

For a long time it was thought that low-energy long- 
wavelength quantum fiuctuations are irrelevant (in d > 1) 
so that Eq. ([4]) is universally trueSl. However, there is 
growing evidence that these 'soft modes' are actually rel- 
evant, changing the form of x or even causing the second- 
order QCP to become weakly first order— i^i^. In this 
paper we restrict ourselves to regimes in which such ef- 
fects is negligible, so that Eq. ^ holds, albeit for some 
v not necessarily equal to 1/2. 

This form for x leads to an effective spin-spin interac- 
tion 



nq 



otq 



(5) 



which diverges at small LOn and q as J Jc] these large 
spin fluctuations are likely to have an important effect 
on superconductivity. 

Eq. ([5]) is a generic expression for the spin-spin interac- 
tion near any kind of ferromagnetic QCP, and should be 
applicable regardless of whether the ferromagnetism and 
superconductivity arise from a single band of itinerant 
electrons or from different bands (or even if the mag- 
netism is due to localized spins, provided that one can 
flnd a model of localized spins that exhibits a second- 
order quantum phase transition) Ultimately, the spin- 
fluctuation-mediated interaction can be represented by a 
low-energy effective action of a form analogous to Eq. 



(6) 



where /j, = / ^^^^ is an integral over the Brillouin zone. 

Note that /^'^(Ck) = IsJ^)^^ where the integral is 
performed over the Fermi surface, weighted by the inverse 
Fermi velocity l/vpCk), where vp{]i.) 



9k 



B. Jacobian T^q for spherical Fermi surface 

The dimensionless effective spin-spin interaction, V^, 
may be calculated in the same way as the effective 
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(^^^ = Q \v + G 

FIG. 2: Diagrams for the matrix self-energy. 



FIG. 1: (Left) 2D circular Fermi surface in k-space. (Right) 
contours of ?)(q) (schematic). 



phonon-mediated attraction in Eq. ([3]): 



J kp / J\z. J Q 



(7) 



where rjq is a Jacobian corresponding to the change of 
variables, which is the 'autocorrelation of the Fermi sur- 
face': 



S{^^) 5(&+q) / / (8) 
k / Jk 



T] can be calculated geometrically in 2D (3D) by con- 
sidering the areas (volumes) of intersection of circular 
annuli (spherical shells). The results are 



m 6(1 



2kF 



'7q 



nqkp \ 1 1 

m e(i 



2k 



' q 

\ '2k p 



2qkp 



(2D), (9) 
(3D). (10) 



Both these expressions are normalized such that rjq = 
V. In both 2D and 3D, it can be shown that jyq goes as 
\/q near q = 0. In both 2D and 3D this, by itself, is an 
integrable singularity. However, when k = and ujn = 0, 
V^q goes as 1/g^ near q = 0. The entire integrand of 
Eq. ([7]) thus goes as ■ We shall show that this makes 
the integral divergent. 



C. Effective interaction near FM-QCP 

For the case of a 2D circular Fermi surface, substituting 
Eq. (O into ^ gives 



. 2 poo 



' r) 



mJKQ 
2Tr^kF Jo 



dq 



Jv 



Kq 



72 

kp 



g2 + + ^ 



K-l- (w„/a)i/3 



(11) 



where ly is the 2D density of states.— Vfj diverges as w^/'^, 
and is cut off for uj„ < ctK^ , at a value of the order of 1/k. 
For a 3D spherical Fermi surface, substituting Eq. dH 
into dZ]) gives 



2kF 



dq 



q- 



aq 



2kf 



2 [kpj ''■«:+(^„/«)i/3 



(12) 



with logarithmic accuracy for small, positive w„. Here j/ 
is the 3D density of states. The integrand has a diver- 
gence cut off by ujn and n^. diverges as In l/w„, and 
is cut off for ujn < ctK^ J at a value of the order of In 1/k. 



D. Strong-coupling equations near FM-QCP 



We now investigate the effect on superconductivity. 
The self-consistent 'strong-coupling' equations for the 
4x4 matrix Matsubara Green functions Q and self- 
energies S are 



k-p P k-p P 



(13) 
(14) 



where p and 7 are suitable generalizations of Pauli 
matrices.25 These equations take into account both 
the scattering and pairing effects of the phonon- 
mediated interaction and the spin-fluctuation- 
mediated-interaction V (although they neglect vertex 
corrections) . The self-energy equation is represented dia- 
grammatically in Fig. [21 For isotropic superconductivity 
with interactions in charge and spin channels the strong- 
coupling equations reduce to the following frequency- 
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where = K ± and C/„, T^, £„, A„ are the 

momentum-integrated ordinary and anomalous Green FIG. 3: Suppression of Tc by static magnetic disorder of 
functions and self-energies. strengtii (dashed line) and by quasi-static spin fluctua- 

tions of strength Vq (solid line). 



E. Abrikosov-Gor'kov equations near FM-QCP 

In the vicinity of the QPT, where k is small, is 
dominated by the divergence at small w„. Now, V!^ is 
sampled at discrete bosonic Matsubara frequencies ujn = 
2TTT{n + 1/2). If the first nonzero Matsubara frequency, 
uji ~ 2ttT, is much larger than the frequency cutoff aK^, 
we can discard all nonzero Matsubara frequencies. This 
corresponds to replacing the dynamic interaction by 
a static one, V^^ ~ Vq Sn-. The spin fluctuations are now 
characterized by a single number, Vg. From Eqs. 
and (HI]), 

Vo^ « J,y (^g^ (^) (2D model), (19) 

« 4^ f ^ V In ^ (3D model). (20) 
2 y kp J K 

Let us also replace the dynamic, phonon-mediated inter- 
action by a BCS-type interaction (V^ = 1/^*-^^), with a 
frequency cutoff lud of the order of the Debyc frequency. 
Then Eqs. US]) and HI]) simplify to 

e„ = En + TV^Gn (21) 

An=T V'^'^^Tn^-TV'^Tn. (22) 

m 

Compare these with the case of static charge and spin 
disorder—, that is, magnetic impurities in a supercon- 
ductor. That system can be described by the BCS model 
with additional random delta-correlated potentials t/^(x) 
and Uf(x) in the charge and spin channels, whose vari- 
ances are and in the sense that (C/'^(x)?7'^(x')) = 
Wdix - x') and (C/f (x);7;(x')) = WSix - x')^y : 

+ 11/BCS/' p^^p^^_ I U'\^)pr^- [ ^[/f(x)5,,, 

(23) 



Let = 2Tri/W'^'^ = where Tc,s are the ordinary 

and spin- flip scattering times, and v is the density of 
states at the Fermi energy. Define = w'^ ± w'^ . One 
then obtains the Abrikosov-Gor'kov equations 

en = en + —Gn (24) 

A„ = T V^^^Tr. + (25) 

Comparing Eqs. ([21]) and with ([M]) and (gS]) shows 
that the quasi-static spin-spin interaction V'' can effec- 
tively be replaced by static spin disorder whose strength 
w'^ = TVq is temperature- dependent. The factor of T 
arises from the Matsubara frequency sum in the Eliash- 
berg equations, and is necessary on dimensional grounds. 

From Abrikosov-Gor'kov theory it is known that 
Tc/TcQ depends only on w'' according to the following 
implicit equation^fi: 

where Tco is the critical temperature of the clean super- 
conductor, and tpiz) is the digamma function. In this 
scenario, Tc falls to zero when = ^ = -^Tco (Fig.[3l 
dashed line) . A finite concentration of magnetic impuri- 
ties is sufficient to destroy superconductivity. However, 
the extra factor of T in the spin-fluctuation scenario in- 
dicates that at lower temperatures there arc fewer exci- 
tations and the depairing effect is less pronounced. Sub- 
stituting in Eq. (j26p . we find 

This is an explicit equation for Tc which is even simpler 
than the AG result (Fig. [31 solid line) ! The limiting forms 
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Tco 

Ik 

Tco 



-V^ for Vo^ ^ 0, 



for Vn' 



oo. 



(28) 
(29) 



Superconductivity is only destroyed when the strength of 
the spin fluctuations becomes inflnitc. Inserting Eq. (HI 
or (gni) into (gl]) gives 



(2D), (30) 
(3D). (31) 



Tc ^ TT 1 
Tr TT 2 



2eT Jv 
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2kF 




K 



In the immediate vicinity of the QCP, the bare in- 
verse correlation length kq is a constant of the order of 
kpi whereas the true inverse correlation length k goes 
to zero according to some power law, k ^ \J — Jd" , 
as J approaches Jc- In Eqs. pO|) and ((3T|) . the dimen- 
sionless quantities Jiy and kp/no are constants of order 
unity, whereas the ratio kq/k (which is the square root of 
the Stoner enhancement factor) becomes infinite. Hence, 
~ \n\j-j^ \ 3D, and |J- J^T in 2D. Note that 
the fluctuations are more severe in 2D than in 3D, and 
Tc is more strongly suppressed. The quantum critical ex- 
ponent I' affects the form of the 2D result but only the 
prefactor of the 3D result. 

It is likely that ( J— Jc) is proportional to (p—pc), where 
p is an experimental control parameter such as pressure. 
Substituting p for J, we obtain predictions that can be 
compared with experiment. 

In the limit k ^ 0, l/ln(l/K) and k are always larger 
than aK^ . Hence our initial assumption 27rT 3> an^ is 
valid in the vicinity of the superconducting transition 
line. This justifles the working above a posteriori. 

Farther away from the QCP, the other factors in 
Eqs. pO|l and ([3T|) may come into play. It is likely that 
changing the experimental parameter p has a significant 
effect not only on J, but also on the degree of nesting of 
the Fermi surface, and hence on kq. Then Tc will deviate 
from a pure logarithm or power law. 
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FIG. 4: Corrections to mean-field phase diagram due to fluc- 
tuations, assuming vertical FM phase boundary. 'N' repre- 
sents the normal, non-magnetic metallic phase. 



eter M attempt to suppress each other. For a single- 
band model in which the same electrons take part in 
Cooper pairing and in 'Stoner' itinerant ferromagnetism, 
the competition is so intense that the coexistence region 
in Fig.|4]is eliminated, leaving a first-order FM/SC phase 
boundary, terminating in a first-order QCP, and the anal- 
ysis in this paper is not relevant. It may be possible, how- 
ever, for FM and SC to coexist in multi-band systems. 
In the following discussion, we presuppose coexistence of 
FM and SC at the mean-field level. 

First let us assume that the magnetic phase bound- 
ary is vertical. On the left of this boundary, M ~ 0, 
so Tc should be constant (dashed line in Fig. 2]). The 
magnetic fluctuations should suppress Tc to zero as de- 
scribed above (solid curve). On the right of the mag- 
netic phase boundary, magnetic order is present, and Tc 
should be suppressed at the mean-fleld level. Suppose 
that the mean-field values of Tc are given by the dashed 
curve in Fig. Now, since the equation for Tc is linear, 
the pair-breaking effect of FM order and FM fluctuations 
combine additively in the argument of the digamma func- 
tion. Since the fluctuations are present on both sides of 
the transition, we expect that Tc{J) is suppressed to zero 
as J ^ Jc from either direction (solid curve) : 



IV. EFFECT OF SDW FLUCTUATIONS ON 



F. Phase diagram 

In the presence of an 'externally applied exchange 
held', such as would arise from artiflcially polarized im- 
purity spins, it is possible for superconducting order to 
exist. In fact, Larkin and Ovchinniko\»ii show that the 
exchange splitting h has no effect on A until it reaches 
about A/2. (The case of an externally applied magnetic 
held is more comphcated due to orbital effects.) 

The case for coexistence of superconductivity and 
spontaneous magnetization is a more difficult one. In a 
fully self-consistent mean-field theory, the superconduct- 
ing order parameter A and the magnetic order param- 



A. r;q near an SDW-QCP 

In order to treat a system close to a SDW instability, 
we must abandon the assumption of a spherical Fermi 
surface. A SDW system typically has a Fermi surface 
which is almost 'nested', meaning that there are por- 
tions of the Fermi surface which almost overlap when 
shifted relative to each other by a nesting vector Q. The 
spin fluctuations with wavevectors close to Q are the low- 
energy quantum critical fluctuations, and it is they that 
play an important role in suppressing superconductivity. 

In order to proceed we need to choose a deflnite model 
for the Fermi surface. Let us begin with a 2D model. 
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FIG. 5: (Left) Circular Fermi surface for < fc^, fcy < ^. 
Solid gridlines indicate fc-space unit cell boundaries; dashed 
lines indicate boundary of 1st Brillouin zone. Empty circle 
is image of Fermi surface under translation by Q. (Right) 
Contours of rjc^ (schematic). 
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FIG. 6: Flattened Fermi surface. The line singularities in 
rjd have a larger magnitude than before. The singularity at 
q = also changes shape, but this is not shown, and is not 
relevant to the discussion. 



Suppose the Fermi surface is a circle of radius kp. Shift- 
ing this by 2kp in any direction gives a circle tangent 
to the original circle, suggesting imperfect nesting. Now 
the anisotropy of the crystal comes into play. Since the 
Fermi surface typically occupies a sizeable fraction of the 
Brillouin zone, we must consider nesting between peri- 
odic images. In k-space, shifting the Fermi surface by 
Q makes it tangent to the original Fermi surface at two 
points rather than one. 

It is useful to perform the analysis in q-space. The 1/q 
singularity in r]q is still present, but there is a more rel- 
evant 1/V2fci? — q singularity on a circle of radius 2k p. 
This circle and its periodic images intersect at the point 
Q (and three other points related to Q by symmetry), 
and there r]q is enhanced by a factor of 2. This construc- 
tion, illustrated in Fig.[5l determines the optimal nesting 
vector—. In summary, the radius of the Fermi surface 
determines the magnitude(s) of the nesting vectors, and 
the lattice anisotropy determines their direction(s). 

Now suppose the Fermi surface undergoes 'pincushion' 
distortion, such that the circular arcs that make up the 
Fermi surface become flatter (Fig. [6]). The optimal vector 
Q may change; more importantly, the 'nesting' singular- 
ities in 77q become larger. Assuming that the Fermi ve- 
locity vp and the density of states v are not significantly 



altered during the distortion. 



'7q 



27rfc, 



2fcc 



' q 

•■2kc 



rnSi^ kc 



where q± is the perpendicular distance from q to the 
Fermi surface in q-space (towards the center of the arc). 
The magnitude of the singularity increases with the 
square root of the radius of the curvature. Remember 
that this is enhanced by a factor of 2 where the singular- 
ities intersect. When kc = oo the Fermi surface becomes 
a square and rj^ cc S{q±). 

In the 3D analogue of this model, the Fermi surface 
consists of eight, approximately triangular, portions of 
spherical surfaces. As the radius of curvature of each 
portion is increased towards infinity, the Fermi surface 
morphs from a sphere to an octahedron: 



2iTkc^Q{l-^) mOiq^) 



Airkf 



Akf 



kc 



(33) 



The discontinuity in r/q is proportional to the radius of 
curvature. 



B. Spin-spin interaction V^q near an SDW-QCP 



We also need an approximation for the spin-spin inter- 
action V* . First consider the bare Green function Xnai 



given by 



/(a) - /(Ck+q) 



(34) 



This integral is somewhat similar to the expression for rjq, 
Eq. (IHI). It is large when q is close to a nesting vector of 
the Fermi surface. In this case, the definition of 'nesting' 
is more stringent: portions of the Fermi surface overlap 
when shifted by the nesting vector Q, and the 'empty' 
side of one portion matches up with the 'filled' side of 
the other. Because of this, Xnq does not have a peak 
near q = 0, unlike ?7q. We shall assume that both Xnq 
and rjq are peaked at the same wavevector Q. This is a 
reasonable assumption if the SDW and SC are caused by 
nesting of the same Fermi surface. 

For small w and for q close to Q, Xnq can be approx- 
imated by a Hertzian of width kq. In this model, it is 
expected that and kq will both depend upon the 
curvature parameter kc- In particular, when kc is large, 
the system is well nested, and kq is expected to be small. 



1/0 2 
XoQ^O 



+ |q - 



(35) 



Then, in X^q, kq is replaced by the 'renormalized' inverse 
spin-fluctuation correlation length k, which is smaller due 
to Stoner enhancement: 



V" = 



(36) 
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FIG. 7: (Left) Electron-hole excitations near the nesting vec- 
tor which give the main contribution to x- (Right) Schematic 
contour plot of xoq (or Voq)- 



where k = kq (1 - JXqqY- Ultimately, the FLEX cal- 
culation requires only the 7i = component, which is a 
Lorentzian in q-space peaked at q = Q: 



V, 



Oq 



|q-QI 



(37) 



C. 



Effective interaction near SDW-QCP; 
suppression of Tc 



The momentum-averaged interaction is given, as be- 
fore, by Eq. ([7|). Since ?7q has a weaker divergence for the 
SDW case than for the FM case, Vq also has a weaker 
divergence, as we shall see below. First consider the 2D 
model. To do the integral over q, we use coordinates 
and q± oriented tangentially and perpendicularly to the 
Fermi surface. There is a factor of 2 arising from the 
overlap of two singularities. 



V = 



(2^ 



k2 +q||2 



X 2 



27r3fcF^ 



dqi 



nl 



, 2 V^C 
TT kp Jo 



/q± J -oo qf + qi_-^ + + ^ 
dq± TT 



8r(5/4)2 



TT ^ 



1/4- 



(38) 



At zero frequency. 



1^ 



(39) 



The strength of the effective interaction increases when 
the Fermi surface becomes flatter or when the coupling 
J is tuned towards Jc- In general, k goes to zero at the 
QCP whereas kc is large but finite. Therefore, near the 
QCP, Vq ^ I/a/kj and by the arguments of the previous 
sections, the superconducting Tc is suppressed as 



Tc 



(40) 



The quantity Jv is of order unity, but kq/ up is likely to 
be small, as explained earlier. Hence the suppression of 
Tc is much weaker for the SDW case than for the FM 
case. 

For the 3D model, an estimate gives 

" J (27r)3 ''1 "1 

1 f^'' f m 
^WrJo dq.Jdq,2nq,^e{q^)kcX 



X 4 



itlJko kc 
kp^ Jo 



dqi 
dq± In 



dq\\ q\\ 



kp 







Vn ~ J I' 1 — — 



kp kp 



2 kp 



(41) 



The integral is not divergent at small q±, so Vq remains 
finite, and Tc is not suppressed to zero. In fact, since 
is small, the suppression should be a small effect. 
Taking into account the terms of order k, we see that 
Vq ^ const — K ~ const — | J — JcT, so Tc ~ const + | J — 
Jcl": the critical temperature has a cusp near the QCP. 

The above arguments apply close to the QCP. Farther 
away from the QCP, the variation of kq may be signifi- 
cant. In fact, it is likely that the quantum phase transi- 
tion is caused by an increase in nesting, so kq decreases 
as the QCP is approached from far. Then the argument 
implies that Tc should increase at first as the QCP is 
approached. This can be justified as follows: although 
Fermi-surface nesting causes an enhancement of x'^ ^J^d 
X, and indirectly enhances the V, it also narrows the 
peak in the interaction in q-space so that the magnetism 
competes less effectively for the Fermi surface. Of course, 
this result is model-dependent. 



V. NUMERICAL RESULTS 

Figure [5] shows numerical results obtained by self- 
consistent solution of Eqs. P^ - (fT5|) . using an Einstein 
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phonon model for V,j and the various approximations de- 
rived above for V*. From the fog-fog pfot (Fig. |8(b)| ) it is 
cfoar that Tc obeys the expected power laws (k and ^/k) 
near a 2D FM or AF transition respectively, and that in 
3D the suppression of Tc is very weak. 



VI. DISCUSSION 

We have made predictions about the analytic behavior 
of the Tc of an isotropic superconductor near ferromag- 
netic or spin-density-wave QCPs: 



2D FM Tc 

2D SDW Tc 

3D FM Tc 

3D SDW Tc 



|J- Jcl" 

\j-JcV'^ 

l/ln(l/K) - l/ln(l/|J- 
const + I J - 



Tc is suppressed all the way to zero by ferromagnetic 
quantum critical fluctuations in 2D and 3D, even before 
magnetic order sets in. Spin-density-wave fluctuations 
have a weaker effect; in 2D they can suppress Tc to zero, 
but in 3D they merely produce a cusp. These predictions 
rely on the assumption of coexistence of magnetism and 
homogeneous, isotropic superconductivity at the mean- 
field level. The results for the 3D FM should be con- 
trasted with the prediction in the p-wave case of a finite 
Tc at the QCPJi 

These results are based on the theory of Abrikosov 
and Gor'kov, which has undergone several refinements 
since its publication in 1961. Calculations including 
higher-order scatterings from a classical spi n^^i^^i^° and 
non-perturbative studies^ suggest that magnetic impu- 
rities lead to states within the superconducting gap; how- 
ever, the predicted Tc (which is all that matters to our 
conclusions) is not qualitatively different from that of 
Abrikosov and Gor'kov. Taking the quantum nature 
of the spins into account^^ suggests that the spins are 
quenched by the Kondo effect at low temperatures (in 
the case of antiferromagnctic exchange J > 0), and thus 
cease to cause pair-breaking, so that superconductivity 
is destroyed only at an infinite concentration of magnetic 
impurities. However, the present paper deals not with 
impurity spins but with quantum critical spin fluctua- 
tions: in a nearly magnetic metal these are not destroyed 
at low temperatures by any Kondo effect, so there is no 
reason to expect such a phenomenon in a nearly magnetic 
superconductor. 

We are not aware of any experimental systems where 
a FM QCP has been seen in a phonon-mediated super- 
conductor. But, as Saxena et al remark in Refi^, very 
few itinerant-electron fcrromagncts studied to date have 



been prepared in a sufficiently pure state, or have been 
'tuned', to be sufficiently close to the border of ferro- 
magnetism, or cooled to sufficiently low temperatures, to 
provide a definitive check of the predictions of theory. 

In Fo22,, signs of superconductivity have been observed 
on the high-pressure side of a first-order structural phase 
transition, but it is unclear if this is (i) conventional 
phonon-mediated superconductivity which is suppressed 
by ferromagnetism in the bcc phase and revealed in the 
hep phase, (ii) unconventional superconductivity medi- 
ated by ferromagnetic spin fluctuations, or (iii) by some 
other type of spin fluctuation present in the hep phase. 
As the transition is flrst-order, the analysis in this paper 
is not applicable. 

The compound MgCNis is a superconductor below at 
8K^, with a BCS-like (s-wave-like) specific heat, and 
there has been a theoretical prediction that it is unstable 
to ferromagnetism upon doping with 12% Na or Li. Such 
a ferromagnetic quantum critical point may be a testing 
ground for this theory. 

In ErRh4B4i and HoMogSg^i^ the FM state de- 
stroys superconductivity at sufficiently low tempera- 
tures, whereas at intermediate temperatures the coex- 
istence gives rise to a compromise oscillatory behavior. 
The physics here is dominated by orbital diamagnetism, 
rather than by paramagnetic spin fluctuations. 

The rare-earth nickel borocarbides^ show coexisting 
SDW and superconductivity, as well as corresponding 
nesting features in the Fermi surface;^^^ Suppression of 
Tc using hydrostatic and chemical pressure to manipu- 
late the exchange interaction has been reported^ but 
the AFM QCP has not been explored. 

In the Bechgaard salt (TMTSF)2PF6, pressure has 
been used to tune the SDW state through a QCP at 
6.4kbar where a superconducting state is found. The 
quasi-2D K-BEDT-TTF salts have a high superconduct- 
ing temperature and are near the boundary of an AF 
insulating phasci^ However, the character and mecha- 
nism of superconductivity in these organic materials is 
still debated, and thus the applicability of our model is 
in doubt. 

To conclude, while there are several candidate model 
systems with which to test these predictions, up to now 
no sufficiently detailed studies have been performed. 
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Strictly speaking, the angular momentum classifications 
s,p,d should be replaced by the appropriate representations 
of the symmetry group of the crystal. 

The Stoner model described above exhibits a quantum crit- 
ical point and a vertical FM phase boundary. The standard 
Heisenberg model, however, gives a finite-temperature 
phase transition whose Curie temperature is proportional 
to J. 

Eq. (Illf) is not a truncated series expansion but rather 
an approximate form with the correct limiting behavior at 
small n and large n. 



